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Abstract 

We demonstrate that Higgs discovery at the LHC is possible in the context of the NMSSM even for 
those scenarios such that the only strongly produced Higgs boson is a very SM-like CP-even scalar 
which decays almost entirely to a pair of relatvely light CP-odd states. In combination with other search 
channels, we are on the verge of demonstrating that detection of at least one of the NMSSM Higgs bosons 
is guaranteed at the LHC for accumulated luminosity of 300 fb^^. 

1. Introduction 

One of the most attractive supersymmetric models is the Next to Minimal Supersymmetric Standard 
Model (NMSSM) (see [1,2] and references therein) which extends the MSSM by the introduction of just 
one singlet superfield, S. When the scalar component of S acquires a TeV scale vacuum expectation 
value (a very natural result in the context of the model), the superpotential term SHuHd generates an 
effective fiHuHd interaction for the Higgs doublet superfields. Such a term is essential for acceptable 
phenomenology. No other SUSY model generates this crucial component of the superpotential in as 
natural a fashion. Thus, the phenomenological implications of the NMSSM at future accelerators should 
be considered very seriously. One aspect of this is the fact that the h, H, A, Higgs sector of the 
MSSM is extended so that there are three CP-even Higgs bosons (/ii,2,3, < nih^ < rrih^), two 
CP-odd Higgs bosons (ai^2> w-ai < "t-^j) (we assume that CP is not violated in the Higgs sector) and 
a charged Higgs pair (h^). An important question is then the extent to which the no-lose theorem for 
MSSM Higgs boson discovery at the LHC (after LEP constraints) is retained when going to the NMSSM; 
i.e. is the LHC guaranteed to find at least one of the /ii,2,3> ai^2> The first exploration of this issue 
appeared in [3], with the conclusion that for substantial portions of parameter space the LHC would be 
unable to detect any of the NMSSM Higgs bosons. Since then, there have been improvements in many of 
the detection modes and the addition of new ones. These will be summarized below and the implications 
reviewed. However, these improvements and additions do not address the possibly important h aa 
type decays that could suppress all other types of signals [3,4]. 

One of the key ingredients in the no-lose theorem for MSSM Higgs boson discovery is the fact 
that relations among the Higgs boson masses are such that decays of the SM-like Higgs boson to AA are 
only possible if ruA is quite small, a region that is ruled out by LEP by virtue of the fact that Z —>■ hA 
pair production was not detected despite the fact that the relevant coupling is large for small m^. In 
the NMSSM, the lighter Higgs bosons, hi or /i2, can be SM-like (in particular being the only Higgs 
with substantial WW/ZZ coupling) without the ai necessarily being heavy. In addition, this situation 
is not excluded by LEP searches for e+e^ —>■ Z* —>■ /ii,2ai since, in the NMSSM, the ai can have 
small Z/i2ai (Zhiai) coupling when hi (/12) is SM-like. [In addition, sum rules require that the Zhiai 
(Zh2ai) coupling is small when the hiWW {h2WW) couplings are near SM strength.] As a result, 
NMSSM parameters that are not excluded by current data can be chosen so that the /ii 2 masses are 
moderate in size (~ 100 — 130 GeV) and the hi —>■ aiai or /12 aici decays are dominant. Dominance 
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of such decays falls outside the scope of the usual detection modes for the SM-like MSSM h on which 
the MSSM no-lose LHC theorem largely relies. 

In Ref. [2], a partial no-lose theorem for NMSSM Higgs boson discovery at the LHC was es- 
tablished. In particular, it was shown that the LHC would be able to detect at least one of the Higgs 
bosons (typically, one of the lighter CP-even Higgs states) throughout the full parameter space of the 
model, excluding only those parameter choices for which there is sensitivity to the model-dependent 
decays of Higgs bosons to other Higgs bosons and/or superparticles. Here, we will address the question 
of whether or not this no-lose theorem can be extended to those regions of NMSSM parameter space for 
which Higgs bosons can decay to other Higgs bosons. We find that the parameter choices such that the 
"standard" discovery modes fail would allow Higgs boson discovery if detection of h ^ aa decays is 
possible. (When used generically, the symbol h will now refer to h = hi, h2 or hs and the symbol a 
will refer to a = oi or 02). Detection of h ^ aa will be difficult since each a will decay primarily to bb 
(or 2 jets if nia < 2m;,), r+r", and, possibly, XiXv yielding final states that will typically have large 
backgrounds at the LHC. 

In [2] we scanned the parameter space, removing parameter choices ruled out by constraints from 
LEP on Higgs boson production, e+e~ Zh or e+e" ha [5], and eliminating parameter choices for 
which one Higgs boson can decay to two other Higgs bosons or a vector boson plus a Higgs boson. For 
the surviving regions of parameter space, we estimated the statistical significances {Nsd = S/y/B) for 
all Higgs boson detection modes so far studied at the LHC [6-9]. These are (with £ = e,iJ,) 

1) 99 h/a 77; 

2) associated Wh/a or tih/a production with 77^=*= in the final state; 

3) associated tt/i/a production with h/a ^ bb; 

4) associated bbh/a production with h/a ^ t'^t~; 

5) gg^h^ ZZ^*^ 4 leptons; 

6) gg-^h~^ WW^*^ e+£~ui?; 

7) WW -^h^ T+T-; 

8) WW -^h^ WW^*l 

For an integrated luminosity of 300 fb~^ at the LHC, all the surviving points yielded Nsd > 10 after 
combining all modes, including the VF-fusion modes. Thus, NMSSM Higgs boson discovery by just one 
detector with L = 300 fb~^ is essentially guaranteed for those portions of parameter space for which 
Higgs boson decays to other Higgs bosons or supersymmetric particles are kinematically forbidden. 

In this work, we investigate the complementary part of the parameter space, where at least one 
Higgs boson decays to other Higgs bosons. To be more precise, we require at least one of the following 
decay modes to be kinematically allowed: 

i) h ^ h'h' , ii) h ^ aa , in) h h^h^ , iv) h ^ aZ , 
v) h h^W^ , vi) a' ^ ha , vii) a ^ hZ , viii) a h^W^ . (1) 

After searching those regions of parameter space for which one or more of the decays i) — viii) is 
allowed, we found that the only subregions for which discovery of a Higgs boson in modes 1) - 8) was 
not possible correspond to NMSSM parameter choices for which (a) there is a light CP-even Higgs boson 
with substantial doublet content that decays mainly to two still lighter CP-odd Higgs states, h aa, and 
(b) all the other Higgs states are either dominantly singlet-like, implying highly suppressed production 
rates, or relatively heavy, decaying to ti, to one of the "difficult" modes i) — viii) or to a pair of sparticles. 
In such cases, the best opportunity for detecting at least one of the NMSSM Higgs bosons is to employ 
WW h production and develop techniques for extracting a signal for the h aa ^ jjT^T~ 
(including jj = bb) process. We have performed a detailed simulation of the aa JJt^t^ final state 
and find that its detection may be possible after accumulating 300 fb~^ in both the ATLAS and CMS 
detectors. 



2, The model and scanning procedures 

We consider the simplest version of the NMSSM [1], where the term 1JLH1H2 in the superpotential of the 
MSSM is replaced by (we use the notation A for the superfield and A for its scalar component field) 

\H^H2S + , (2) 

so that the superpotential is scale invariant. We make no assumption on "universal" soft terms. Hence, 
the five soft supersymmetry breaking terms 

m\^Hl + mji^Hl + m^S^ + \AxH1H2S + (3) 

are considered as independent. The masses and/or couplings of sparticles will be such that their contri- 
butions to the loop diagrams inducing Higgs boson production by gluon fusion and Higgs boson decay 
into 77 are negligible. In the gaugino sector, we chose M2 = 1 TeV (at low scales). Assuming universal 
gaugino masses at the coupling constant unification scale, this yields Mi ^ 500 GeV and M3 ^ 3 TcV. 
In the squark sector, as particularly relevant for the top squarks which appear in the radiative corrections 
to the Higgs potential, we chose the soft masses mq = rriT = Mgusy = 1 TeV, and varied the stop 
mixing parameter 

X = 2 ^ II — 

As in the MSSM, the value = \/6 - so called maximal mixing - maximizes the radiative corrections 
to the Higgs boson masses, and we found that it leads to the most challenging points in the parameter 
space of the NMSSM. We adopt the convention X,k>0, in which tan (3 can have either sign. We require 
l/Lieffl > 100 GeV; otherwise a light chargino would have been detected at LEP. The only possibly light 
SUSY particle will be the x?. A Ught Xi is a frequent characteristic of parameter choices that yield a 
Ught ai. 

We have performed a numerical scan over the free parameters. For each point, we computed 
the masses and mixings of the CP-even and CP-odd Higgs bosons, hi {i = 1,2, 3) and aj (J = 1, 2), 
taking into account radiative corrections up to the dominant two loop terms, as described in [10]. We 
eUminated parameter choices excluded by LEP constraints [5] on e~^e~ Zhi and e+e~ — ^ hiaj. The 
latter provides an upper bound on the ZhiUj reduced coupling, R'-, as a function of rrih^ + ruaj for 
m/j. ~ niaj. Finally, we calculated m/j± and required m/j± > 155 GeV, so that t ^ h^h would not be 
seen. 

In order to probe the complementary part of the parameter space as compared to the scanning of 
Ref. [2], we required that at least one of the decay modes i) — viii) is allowed. For each Higgs state, 
we calculated all branching ratios including those for modes i) — viii), using an adapted version of the 
FORTRAN code HDECAY [11]. We then estimated the expected statistical significances at the LHC in 
all Higgs boson detection modes 1) - 8) by rescaling results for the SM Higgs boson and/or the MSSM 
h, H and/or A. The rescaling factors are determined by Ri, ti and bi = ti, the ratios of the VVhi, tthi 
and bbhi,T^T~hi couplings, respectively, to those of a SM Higgs boson. Of course < 1, but U 
and bi can be larger, smaller or even differ in sign with respect to the SM. For the CP-odd Higgs bosons, 
i?^ = at tree-level; t'j and bj are the ratios of the ij^ coupfings for tt and bb, respectively, relative 
to SM-like strength. A detailed discussion of the procedures for rescaling SM and MSSM simulation 
results for the statistical significances in channels 1) - 8) will appear elsewhere. 

In our set of randomly scanned points, we selected those for which all the statistical significances 
in modes 1) - 8) are below 5a. We obtained a lot of points, all with similar characteristics. Namely, in 
the Higgs spectrum, we always have a very SM-Uke CP-even Higgs boson with a mass between 115 and 
135 GeV (i.e. above the LEP limit), which can be either hi or /12, with a reduced coupling to the gauge 
bosons i?i ~ 1 or i?2 — 1, respectively. This state decays dominantly to a pair of (very) Ught CP-odd 
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Table 1: Properties of selected scenarios that could escape detection at the LHC. In the table, Ri — QhiVV / QHsmVY , ti = 
Shitt/dhsMti = Sh-bb/ShsMbb '^hsM = rrih-; t'l and b'l are the 175 couplings of ai to ti and bb normalized 

relative to the scalar ti and bb SM Higgs couplings. We also give the gg fusion production rate ratio, gg — > hi/gg — » hsM, 
for ruhsM = • Important absolute branching ratios are displayed. For points 2 and 6, the decays oi — > jj (j ^ b) 
have BR{ai jj) ~ 1 — BR{ai r^r~). For the heavy hs and az, we give only their masses. For all points 1-6, 
the statistical significances for the detection of any Higgs boson in any of the channels 1) - 8) are tiny; the next-to-last row 
gives their maximum together with the process number and the corresponding Higgs state. The last row gives the statistical 
significance of the new WW h ^ aa ^ jjT~^T~ [h = hi {h = /12) for points 1-3 (4-6)] LHC signal explored here. 



states, oioi, with ttIq^ between 5 and 65 GeV. The singlet component of ai cannot be dominant if we are 
to have a large hi —>■ oiai or /12 —>■ aiai branching ratio when the hi or /12, respectively, is the SM-like 
Higgs boson. Further, when the hi or /12 is very SM-like, one has small Zhiai or Zh2ai coupling, 
respectively, so that e+e" — > hiai or e"*"e~ /i2ai associated production places no constraint on the 
light CP-odd state at LEP. We have selected six difficult benchmark points, displayed in Tabled These 
are such that oi XiXi decays are negligible or forbidden. (Techniques for cases such that x\x\ decay 
modes are important are under development.) For points 1 - 3, ^1 is the SM-like CP-even state, while 
for points 4 - 6 it is /i2. We have selected the points so that there is some variation in the /ii 2 and 
ai masses. The main characteristics of the benchmark points are displayed in Tabled Note the large 
BR{h —>■ aiai) of the SM-like h (h = hi for points 1-3 and h = h2 for points 4 -6). For points 4-6, 
with rrifi^ < 100 GeV, the hi is mainly singlet. As a result, the Zhiai coupling is very small, implying 
no LEP constraints on the hi and ai from e^e^ — > hiai production. 

We note that in the case of the points 1-3, the /i2 would not be detectable either at the LHC or 
at a Linear Collider (LC). For points 4-6, the hi, though light, is singlet in nature and would not be 
detectable. Further, the h^ or 02 will only be detectable for points 1 - 6 if a super high energy LC is 
eventually built so that e+e~ Z ^ h^a2 is possible. Thus, we will focus on searching for the SM-like 
hi (/i2) for points 1 - 3 (4 - 6) using the dominant hi (/i2) —>■ aiai decay mode. 

In the case of points 2 and 6, the ai t^t^ decays are dominant. The final state of interest will 
be jjr^T^, where the jj actually comes primarily from aiai — > t~^t~t~^t~ followed by jet decays of 
two of the r's: r+r^ jj + u's. (The contribution from direct ai jj decays to the jjr^r^ final 
state is relatively small for points 2 and 6.) In what follows, when we speak of r+r^, we refer to those 
r's that are seen in the r+r^ —>■ £^£^ + v's final state (£ = e, /i). For points 1 and 3-5 BR{ai —>■ bb) 
is substantial. The relevant final state is bbT^T~ . Nonetheless, we begin with a study of the backgrounds 
and signals without requiring 6-tagging. With our latest cuts, we will see that 6-tagging is not necessary 
to overcome the apriori large Drell-Yan r+r^-f-jets background. It is eliminated by stringent cuts for 
finding the highly energetic forward / backward jets characteristic of the WW the fusion process. As a 
result, we will find good signals for all 6 of our points. 

In principle, one could explore final states other than bbr^r^ (or jjr^r^ for points 2 and 6). 
However, all other channels will be much more problematical at the LHC. A 46-signal would be bur- 
dened by a large QCD background even after implementing ^-tagging. A 4j-signal would be completely 
swamped by QCD background. Meanwhile, the 4r-channel (by which we mean that all taus decay 
leptonically) would not allow one to reconstruct the hi, ^2 resonances. 

In the case of the 262r (or 2j2T) signature, we identify the r's through their leptonic decays to 
electrons and muons. Thus, they will yield some amount of missing (transverse) momentum, p^^gg. 
This missing transverse momentum can be projected onto the visible e, //-momenta in an attempt to 
reconstruct the parent r-direction. 

3. Monte Carlo Results for the LHC 

Let us now focus on the WW —^h—>^aa channel that we believe provides the best hope for Higgs 
detection in these difficult NMSSM cases. (We reemphasize that the hi [cases 1 - 3] or /i2 [cases 4 - 
6] has nearly full SM strength coupling to WW.) The bbr^r^ (or 2jr+r^, for points 2 and 6) final 
state of relevance is complex and subject to large backgrounds, and the ai masses of interest are very 
modest in size. In order to extract the WW fusion 2j2r NMSSM Higgs boson signature, it is crucial 
to strongly exploit forward and backward jet tagging on the light quarks emerging after the double W- 
strahlung preceding Typi^-fusion. We also require two additional central jets (from one of the a's) and 
two opposite sign central leptons (£ = e, /x) coming from the the t^t~ emerging from the decay of the 
other a. By imposing stringent forward / backward jet tagging cuts, we remove the otherwise very large 
background from Drell-Yan r+r^ + jets production. In the end, the most important background is due 
to ti production and decay via the purely SM process, gg ^ ti —>■ bbW~^W~ —>■ bbr^T^ + p^i^s' 



association with forward and backward jet radiation. 

We have employed numerical simulations based on a version of HERWIG v6 . 4 [12-14] modified 
to allow for appropriate NMSSM couplings and decay rates. Calorimeter emulation was performed using 
the GET JET code [15]. Since the ai will not have been detected previously, we must assume a value 
for niai . In dealing with actual experimental data, it will be necessary to repeat the analysis for densely 
spaced m^i values and look for the choice that produces the best signal. We look among the central 
jets for the combination with invariant mass Mjj closest to rua^. In Fig. ^ we show the Mjj^+^- 
invariant mass distribution obtained after cuts, but before 6-tagging or inclusion of K factors — the plot 
presented assumes that we have hit on the correct rua^ choice. 
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Fig. 1: We plot da / dMjj^+^- [fb/10 GeV] vs Afjj,+ ,- [GeV] for signals and backgrounds after basic event selections, but 
before h tagging. The lines corresponding to points 4 and 5 are visually indistinguishable. No K factors are included. 

The selection strategy adopted is a more refined (as regards forward / backward jet tagging) version 
of that summarized in [16]. It is clearly efficient in reconstructing the hi (for points 1-3) and /12 (for 
points 4-6) masses from the jjT~^T~ system, as one can appreciate by noting the peaks appearing at 
Mjj^+^- « 100 GeV. In contrast, the heavy Higgs resonances at for points 1-3 and the rather light 
resonances at m^^ for points 4-6 (recall Table Q do not appear, the former mainly because of the very 
poor production rates and the latter due to the fact that either the hi — > oiai decay mode is not open 
(points 4, 5) or - if it is - the jets and e//i-leptons eventually emerging from the ai decays are too soft 
to pass the acceptance cuts (point 6, for which = 7 GeV and m/jj = 51 GeV). For all six NMSSM 
setups, the Higgs resonance produces a bump below the end of the low mass tail of the tt background 
(see the insert in Fig.Q. Note how small the DY r+r^ background is after strong forward / backward 
jet tagging. Since the main surviving background is from tt production, b tagging is not helpful. For 
points 2 and 6, for which the signal has no 5's in the final state, anti- 6- tagging might be useful, but has 
not been considered here. 

To estimate S/y/B, we assume L = 300 fb~\ a factor of 1.1 for the WW fusion signal and 
K factors of 1, 1 and 1.6 for the DY r+T~, ZZ production and tt backgrounds, respectively. (These 
K factors are not included in the plot of Fig. [l]) We sum events over the region 40 < M, ^y+t- ^ 




150 GeV. (Had we only included masses below 130 GeV, we would have had no ti background, and the 
S/^/B values would be enormous. However, we are concerned that this absence of tt background below 
130 GeV might be a reflection of limited Monte Carlo statistics. As a result we have taken the more 
conservative approach of at least including the first few bins for which our Monte Carlo does predict 
some ti background.) 

For points 1, 2, 3, 4, 5 and 6, we obtain signal rates of about S = 1636, 702, 2235, 2041, 2013, 
and 683, respectively. The tt+jets background rate is Bu ~ 795. The ZZ background rate is Bzz ~ 6. 
The DY T^T~ background rate is negligible. (We are continuing to increase our statistics to get a fully 
reliable estimate.) The resulting Ngo = S/\/~B values for points 1-6 are 50, 22, 69, 63, 62, and 21, 
respectively. The smaller values for points 2 and 6 are simply a reflection of the difficulty of isolating 
and reconstructing the two jets coming from the decay of a very Ught ai. Overall, these preUminary 
results are very encouraging and suggest that a no-lose theorem for NMSSM Higgs detection at the LHC 
is close at hand. 

4. Conclusions 

In summary, we have obtained a statistically very significant LHC signal in the jjr^T^ final state of 
WW fusion for cases in which the NMSSM parameters are such that the most SM-like of the CP-even 
Higgs bosons, h, is relatively light and decays primarily to a pair of CP-odd Higgs states, h ^ aa with 
a bb,T~^T~ if iria > 2mi, or a ^ H^t'^t' if rua < 2mb. The statistical significances are (at least) 
of order 50 to 70 for points with rUa > 2m,b and of order 20 for points with < 2mb. These high 
significances were obtained by imposing stringent cuts requiring highly energetic forward / backward jets 
in order to isolate the WW fusion signal process from backgrounds such as DY t'^t^ pair production. 
Still, this signal will be the only evidence for Higgs bosons at the LHC. A future LC will probably be 
essential in order to confirm that the enhancement seen at the LHC really does correspond to a Higgs 
boson. At the LC, discovery of a light SM-like h is guaranteed to be possible in the Zh final state using 
the recoil mass technique [17]. 

In the present study, we have not explored the cases in which the oi XiXi decay has a large 
branching ratio. Detecting a Higgs signal in such cases will require a rather different procedure. Work 

on the WW h ^ invisible signal is in progress [18]. 

As we have stressed, for parameter space points of the type we have discussed here, detection of 
any of the other MSSM Higgs bosons is likely to be impossible at the LHC and is likely to require an LC 
with ^Sg+g- above the relevant thresholds for h'a' production, where h' and a' are heavy CP-even and 
CP-odd Higgs bosons, respectively. 

Although results for the LHC indicate that Higgs boson discovery will be possible for the type of 
situations we have considered, it is clearly important to refine and improve the techniques for extracting 
a signal. This will almost certainly be possible once data is in hand and the tt background can be more 
completely modeled. 

Clearly, if SUSY is discovered and WW WW scattering is found to be perturbative at WW 
energies of 1 TeV (and higher), and yet no Higgs bosons are detected in the standard MSSM modes, a 
careful search for the signal we have considered should have a high priority. 

Finally, we should remark that the h ^ aa search channel considered here in the NMSSM frame- 
work is also highly relevant for a general two-Higgs-doublet model, 2HDM. It is really quite possible 
that the most SM-like CP-even Higgs boson of a 2HDM will decay primarily to two CP-odd states. This 
is possible even if the CP-even state is quite heavy, unlike the NMSSM cases considered here. If CP 
violation is introduced in the Higgs sector, either at tree-level or as a result of one-loop corrections (as, 
for example, is possible in the MSSM), h h'h" decays will generally be present. The critical signal 
will be the same as that considered here. 
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